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A universal differential equation is a nontriviai differential equation the solutions of which ap- 
proximate to arbitrary accuracy any continuous function on any interval of the real line. On the 
other hand, there has been much interest in exactly solvable chaotic maps. An important problem 
is to generalize these results to continuous systems. Theoretical analysis would allow us to prove 
theorems about these systems and predict new phenomena. In the present paper we discuss the 
concept of universal functions and their relevance to the theory of universal differential equations. 
We present a connection between universal functions and solutions to chaotic systems. We will 
show the statistical independence between X{t) and X{t + t) (when r is not equal to zero) and 
X{t) is a solution to some chaotic systems. We will construct universal functions that behave as 
^ ■ delta-correlated noise. We will construct universal dynamical systems with truly noisy solutions. 

O ' We will discuss physically realizable dynamical systems with universal-like properties. 
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I. INTRODUCTION 



Recently there has been great interest in exactly solvable chaotic systems [E S S 0, H, S 0|- ^- Ulam and J. 
von Neumann were the first to prove that the general solution to the logistic map can be found [J, It is very 
^ I important to extend these results to continuous systems. Theoretical analysis would allow us to prove theorems about 
these systems and predict new phenomena. Another surprising parallel development is that of "universal differential 
equations" . A universal differential equation is a nontriviai differential-algebraic equation with the property that 
its solutions a ppr oximate to arbitrary accuracy any continuous function on any interval of the real line [1, [lol . [Til , 
m, [13, [3j [TsI. I16|| . Rubel found the first known universal differential equation by showing that there are differential 
^ I equations of low order (e.g. fourth-order) which have solutions arbitrarily close to any prescribed function The 
existence of universal differential equations illustrate the amazing complexity that solutions of low-order dynamical 
systems can have. In the present paper, we will review some developments in the areas of universal differential 
equations and exactly solvable chaotic dynamical systems. We will discuss the concept of "universal functions" and 
their relevance to the theory of universal differential equations. We will show the statistical independence between 
X{t) and X{t + t) (when r 7^ 0) and X{t) is the solution to some chaotic dynamical systems. We will present a 
connection between universal functions and solutions to chaotic systems. We will construct universal functions that 
, behave as ^-correlated noise. We will construct universal differential equations with truly noisy solutions. We will 
discuss physically realizable dynamical systems with universal properties and their potential applications in secure 
communications and analog computing. 

>< 

; II. UNIVERSAL DIFFERENTIAL EQUATIONS 

Rubel's theorem [1] is: 

There exists a nontriviai fourth- order differential equation 

P{y',y",y"',y"")^0, (1) 

where y = P is polynomial in four variables, with integer coefficients, such that for any continuous function (f>{t) 
on (— cxd,cxd) and for any positive continuous function s{t) on (— cxd,oo), there exists a C°° solution y(t) such that 
\y{t) — (/'(t)! < s{t) for all t on (— cx3,oo). A particular example of Eq.ijT]) is the following: 

3 y''y'Y"' - 4 y'^'^y"" + 6 y''y"'y"'y"" + 2Ay''y"'y"" 
- 29/y"'2/"'' - I2y''y"y"'' + I2y"' = 0, 
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FIG. 1: Solution y{t) for (a) equatfon ([3]), and (b) equatfon ^ with to = 4. 



DufRn ^IQi] has found two additional families of universal differential equations: 



my'^y"" + (2 - 3m)y'y"y"' + 2(m - l)y"^ = 



and 



m^y'^y"" + 3m(l - ra)y'y"y"' + {2m^ - 3m + l)y"'' = 



where to > 3. 

Recently, Briggs has found a new family: 



y"y"" - ?,y'y"y"' + 2(1 - n-^)y"' ^ 0, 



(3) 
(4) 

(5) 



for n > 3. 

The solutions are C". 

We would like to make some observations about these equations. 

Rubel's function y{t) in Eq.Q is C°° but not real-analytic, typically having a countable number of essential 
singularities. 

The functions used to reconstruct the solutions to Eq.([2) are of the form 



y Af{at + P) + B, 
for -1 < t < 1, with g{t) ^ for all other t. 



(6) 



where f{t)^ J g{t)dt, g{t) = exp 

The solutions to equations ^ are trigonometric splines. The kernel g{t) is defined as 

= a[cos6t + c]™. (7) 

On the other hand, the kernel employed to obtain the solutions to Eq. (j4]) are polynomial splines: 

g{t) = a[l-{ht + cfr. (8) 



The solutions are very unstable. This phenomenon can be observed in Fig.(|Ta| and Fig. ([lb 



III. EXACTLY SOLVABLE CHAOTIC MAPS 

Ulam and von Neumann T, ^ proved that the function Xn = sin^(6l7r2") is the general solution to the logistic map 

X„+i = 4X„(1 - X„). (9) 

Recently, many papers have been dedicated to exactly solvable chaotic maps [1, 0, d, H, 0, El, El, [H, [2l|. In 
some of these papers, the authors not only find the explicit functions X„ that solve the maps, but also they discuss 
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the statistical properties of the sequences generated by the chaotic maps. For many of the maps discussed in the 
mentioned papers, the exact solution can be written as 

X„ = P(0fc"), (10) 

where P{t) is a periodic function, is a fixed real parameter and k is an integer. 

For instance, Xn = cos(27r0 3") is the solution to map Xn+i = X„(4X„^ — 3). This is a particular case of the 
Chebyshev maps 0. 

Very interesting chaotic maps of type Xn+i = /(X„) can be constructed when function P(t) is a combination of 
trigonometric, elliptic and other generalized periodic functions 



IV. STATISTICAL INDEPENDENCE 



We will consider statistical independence in the sense of M. Kac |22l |23|. \24. l2a|. In this framework, two functions 
.fi{t), .f2{t) are independent if the proportion of time during which, simultaneously, fi{t) < ai and f2{i) < ^2 is equal 
to the product of the proportions of time during which separately fi{t) < ai, f2{t) < ct2- 

First, we will present some results about functions of natural argument n. 

Consider a generalization to the functions that are exact solutions to the Chebyshev maps: 

= cos(27r6lZ"), (11) 

where Z is a generic real number. 

Any set of subsequences Xs,Xs^i,--- ,Xs^r (for any r) constitutes a set of statistically independent random 
variables. 

li E{X) is the expected value of quantity X, let us define the r- order correlations [18.] : 

1 

E{Xn,Xn, ■ • •X„J - / dXa[p{Xa)Xn,Xn, ■ • -X^J. (12) 



-1 

Here i;(X„) = 0, -1 ^ X„ «C 1, p(X) = = cos(2^0). 

In Ref. [13] it is shown that 

E{x-^x:i, ■ ■ ■ x:i,) ^ E{xT)E{x:i,) ■ ■ ■ e^xx) (13) 

for all positive integers Uq, ni, 712, • • • ,71^. 

The results about the independence of subsequences of function X„ — cos{2tt9Z") can be extended to more general 
functions as the following: 

x„-p(0rz"), (14) 



where P{t) is a periodic function, T is the period of P{t) and 6* is a parameter |27|. 

M. Kac [13, m, [13i HE] has studied independence of different continuous functions, e.g. fi{t) = cos{t),f2{t) = 
cos{\/2t). However, these functions are periodic [2^ [23l. [2^. [25l. [26| . 

We are interested in the independence of continuous functions in the sense that the same function can produce 
statistically independent random variables if evaluated at different times. That is, the functions /(<) and /(< + t) 
should be statistically independent if r 0. 

Periodic functions will never possess this property. 

Using the theorems of Refs. [1^ [1H|, we obtain the result that functions /i = cos(Ait), /2 = cos(A2i), ■ • ■ , fr = 
cos(Ari) constitute a set of independent functions when the numbers Ai, A2, • • • , Ar are linearly independent over the 
rationals. This is equivalent to the fact that if ai,a2, ■ ■■ ,ar are rational, aiAi + a2A2 + • • • + a-rAr = only if 
tti, a2, • ■ ■ yCtr are all zero. 

Consider the following continuous function constructed as a continuous analogous to the function (11): 

X{t) = cosiOe''*). (15) 

It is not difficult to see that the functions X{t) and Y{t) = X{t + r) are independent (in M. Kac's sense) for 
(Lebesgue) almost all r > 0. Note that Y{s) = cos As, X{s) ~ coss, where s — fle***, A = e*""^. 

T he p roof of the fact that Y{s) = cos(As) and X{s) = coss are independent is trivial based on the theorems of M. 
Kac[ll[ll,[li,[13. 
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V. ALGEBRAIC DIFFERENTIAL EQUATIONS AND CHAOTIC SOLUTIONS 



There exists a nontrivial fourth- order differential equations of type: 

Piy',y",y"',y"")^0, (16) 

where P is a polynomial in four variables, with integer coefficients, such that there exist solutions y{t) with the 
property that y(t) and y(t + r) are statistically independent functions for most r. 
A particular example of Eq. p6)) is the following: 

2y'y" - 3y'y"' + y'y"" - y"y"' + {y"f = 0. (17) 
The function y(t) = cos(e^) is a solution to this equation. 

Here we should make a remark. The differential equations 121), ([3]), ([1]), ((S]), and (fTT]) are similar in the sense that 
they are all of the type ([T]) or p6|) . and that all the terms are nonlinear. 

We have shown that the equation ([T7]) possesses entire solutions as the following y{t) = cos(e*) which has the 
property that y{t) and ?/(t + T) are not correlated. This solution is real-analytic on (—00,00). This result shows that 
equations of type (I16p can generate high complexity. 



VI. UNIVERSAL FUNCTIONS 

In this section we will follow the approach to universal equations developed by R. C. Buck Q and M. Boshernitzan 
[1^. In this paper, a family of functions will be called universal if these functions are dense in the space C[i] of 
all continuous real functions on the interval /C M. R. C. Buck has obtained universal partial algebraic differential 
equations using a very deep method: Kolmogorov's solution of Hilbert's Thirteen Problem. He has found that solutions 
to a smooth PDE can be dense in C[i][9]. 

The advantage of this approach is that we can construct universal functions which are real- analytic on R — (—00, 00). 
Thus, the corresponding universal equations will have solutions that are real- analytic on R = (—00, 00). 

For any interval Jc R = (—00,00), C[i] denotes the Banach space of real- valued continuous bounded functions on 

/. 

Boshernitzan has studied the following families of functions: 

t+a 

/bd 
—— cos{e')ds + c, (18) 

1 + — cos(os) 



where d > 0, a, 6, and c are real parameters. 

Another important family of functions is defined as follows: 

t+a 

y{t)=bn J cos2"'(6s)cos(e'*)ds-Hc, (19) 


where a, b, c are parameters and n > 1 is an integer constant. 

We should stress that all the functions in the family given by Eq. (|19p are real-analytic on R = (—00, 00) and entire. 

Boshernitzan has proved that each of the families of functions, (|18p and p9| . is dense in C[I\, for any compact 
interval /= [a, b] C R. 

So these functions are universal. Hence it is possible to construct universal systems using these universal families 
of functions. 

Other results of Boshernitzan are the following. There exists a nontrivial sixth- order algebraic differential equation 
of the form 

P{y',y",--- ,y^'^)=0, (20) 

such that any functions in the family (fT5)) is a solution. And there exists a nontrivial seventh- order algebraic differential 
equation of the form 

P{y',y",--- ,y^'^) = 0, (21) 
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FIG. 2: Time- series generated by ((T9| . 



such that any function in the family (jl9[) is a solution. 

From these theorems one can obtain an important results for us: There exists a nontrivial seventh- order differential 
equation, the real- analytic entire solutions of which are dense in C[i], for any compact interval Jc M. Thus this 
equation is universal. Boshernitzan has not constructed this equation explicitly. 

The possibility of an entire approximation for any dynamics is very promising for practical applications. 



VII. UNIVERSAL SYSTEMS OF DIFFERENTIAL EQUATIONS 

Using the inverse problem techniques and theorems of the works H, [13, [ll|, [l^, [H, H, [iM [3i ^Lnd the results 
contained in the present paper, it is possible to write down a system of differential equations such that any function 
of the family (I19p is a solution: 

Pi[x,x',--- ,x"']^Q., (22) 

P2[y.y'r-- ,y"']^Q, (23) 

z' = Axy. (24) 

Eq. (|22|) is constructed in such a way that all the functions x = a[cos(6a;-f c)]'" are solutions. The specific polynomial 
Pi is 

Pi = mx" x? — (3m — 2)x"x'x H- (2m — 2)x'^ ^ (25) 

where m = 2n? > 2. 

Eq. ((23|) is constructed in such a way that y{t) = cos(e*) is a solution. The specific polynomial P2 is defined as 

P2 = yy'" - iyy" - yy" + y"^ + 2yy. (26) 

The vector solution (x, z) to the system of equations (|22II24I) is such that for the variable z(t) the functions ([TO]) are 
solutions. 

Thus, the system of differential equations ()22ll24p can generate, in variable z(i), universal functions. 
Note that this is a seventh- order dynamical system as was predicted by Boshernitzan [T3| . 
The system of equations (|22j|24p is presented in this paper for the first time. 



VIII. NOISY FUNCTIONS 



In this section we use several conce pts and results from probability theory and mathematical statistics which can 
be consulted for instance in the books [28l [29| . 

In sections (jIVp and (jV| we discussed the function y[t) — cos(e*). Note that the statistical independence between 
functions yi(t) = cos(e*) and y2{t) = cos(e*+'^) implies the following relationship: 



E[y'l^{t)ytit)]^E[y'[^it)]E[y!^^it)], 



(27) 
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for all positive integers fci and fc2. 

Here E[x{t)] is the expected value of quantity x{t). It can be calculated as follows: 

T 

E[xit)] = lim i / x{t)dt. (28) 



As E[yi{t)] — 0, we obtain that 

E[y{t)yit')] = 0, (29) 

for t ^ t'. 

In fact, a direct calculation of the autocorrelation function confirms this results: 

T 

C(t) = lim - / cos(e*)cos(e*+^)dt = (30) 
r-.oo T J 



for T^O. 

Another important statistical property of the independent functions yi{t) and y2{t) is 



where r]{y) is the probability density. 
In fact, 



v{yi{t),y2{t)) = v{yi{t))v{y2{t)), (3i) 



77(2/1) - i — ^ , (32) 



ry(y2) = ^ , (33) 



v{yi,y2)^ . ==■ (34) 



Fig. ([2]) shows these properties. 

Let us introduce a generalized function 



yit) = cos(0(i)). (35) 



We should remark here that the argument of function ([SS)) . (pit), does not need to be exponential all the time for 
t ^ 00, in order to generate noisy dynamics. 

In fact, it is sufficient for function (j)(t) to be a bounded nonperiodic oscillating function which possesses repeating 
intervals with truncated exponential behavior's^. 

A deep analysis of the Boshernitzan's proofs of the fact that the families functions (I18|) and (fT9|) are dense in C[l] 
shows that the random behavior of functions of type ([55]) is crucial [Tsj. 

We are going to present here two examples of these functions. The first is defined by the equation 

x{t) ~ cos{Aexp [ai sin(a;it + (j>i) + 02 sin(tJ2i + ^2) + 03 sin(a;3t + (ps)]} . (36) 

The second function is given by the equation 

y{t) = cos sinh [ai cos(ijJit + (pi) +02 cos(a;2i + ^2)] 
+B2 cosh [fla cos{uj-it + (p-^) + 04 cos(ijJ4t + ^4)]} . 

Note that if the frequencies uJi are linearly independent over the rationals, in both cases, the argument function 
(p{t) is a nonperiodic function with truncated exponential behavior. 

Theoretical and numerical investigatio ns g ive the result that, when the parameters satisfy certain conditions, they 
behave as the solutions to chaotic systems [20!] . Figure (H]) shows that the time- series generated by is very complex. 
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FIG. 4: Time- series generated by function ([3^ 



Furthermore, they behave in such a way that yi = y{t) and 7/2 = y{t + T) are statistically independent functions (in 
M. Kac's sense) for t 7^ 0. An illustration of these properties is that 



Other noisy functions can be obtained using another generalization 

x{t) = p[m, 



(38) 



(39) 
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where P{y) is a general periodic function, and (j){t) is a nonperiodic exponential-like function as before. 
The probability density of these functions depends on the choice of P{y)- 
An important example is the following 

a;(t) ln{tan2[0(t)]} . (40) 

The probability density of the time- series produced by function (|40|) is a Gaussian-like law[30l|. 
Another remarkable property of function (I40p is the following: 

E[x{t)x{t')] = D6{t-t'), (41) 

where E[x{t)x{t')] is defined as in Eq. (|27|) . and d{t — t') is Dirac's delta-function. 

Thus, function (HDl) possesses all the properties usually required in the stochastic equations with 6- correlated noisy 
perturbations (31I |32| . 



IX. PROPERTIES OF CHAOTIC SOLUTIONS 



A Lyapunov exponent of a dynamical system characterizes the rate of separation of infinitesimally close trajectories. 
Suppose Szq is the initial distance between two trajectories, and 6z{t) is the distance between the trajectories at 
time t. 

The maximal Lyapunov exponent can be defined as follows: 



A = lim — In 

t — >oc t 



Sz{t) 



5zo 



(42) 



Consider the dynamical system 



The variational equation is 



§-Fix). (43) 



^ - i?.F(f)0(i). (44) 



The Lyapunov exponent A satisfies the equation 



A ^ lim -ln|0(i,a;o)| . (45) 

t^oo t 

A solution will be considered chaotic if A > 0. 

A related property of chaotic systems is sensitive dependence on initial conditions [ssl [s^l . 
Function 

y{t) = cos[exp(i -I- (/>)] (46) 

has been shown to be the solution to some dynamical system. Here (j) can define the initial condition. 
The Lyapunov exponent of solution can be calculated exactly A = In e = 1 > 0. 

The functions ([M)) . ([57)) . and (UOl) are also chaotic solutions in this sense. The initial conditions are defined by the 
parameters ipi, 02 and (^3. 

Let us discuss sensitive dependence on initial conditions in the context of these functions. 

Let S be the set of functions defined by one of the families given by equations ([36]) . (f37|) . (I40|) . and ((46|) . 

A set S exhibits sensitive dependence if there is a S such that for any e > and each 2/1(0, t) in S, there is a y2{4>\ t), 
also in S, such that \yi{(f>,0) — 2/2(0', 0)| < e, and \yi{(f>,t) — j/2(0',i)| > S for some t > 0. 

The exponential behavior in the arguments of these functions ( (|36|) . (f37| . (j40|) . and ((46|l ) makes them chaotic (see 
a discussion in Ref. [s^). 

All these functions possess equivalent dynamical and statistical properties. 

Following Boshernitzan theory [l^ of modulo 1 sequences, it is possible to construct the following families of 
universal functions: 

t+a 

y{t) = bn I cos^''\bs)x{s)ds + c, (47) 



where x{s) is one of the functions ( ([M)) . ([57]) . (HDl) ) . 
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X. DYNAMICAL SYSTEMS WITH NOISY SOLUTIONS 

In this section we will construct autonomous dynamical systems, the solutions of which are the noisy functions 
discussed in the previous section. 

Consider the following dynamical system 



x[ — Xi[l — {xi - 


^2/1)] - wiyi, 


(48) 


y'l = 2/1 [1 - (xl H 


h yl)] + ujixi, 


(49) 


4 = X2[i - (4 - 


^2/2)] - W22/2, 


(50) 


2/2 = 2/2[l - (4 H 




(51) 


4 = [1 - (4 - 


^2/3)] - W32/3, 


(52) 


2/3 ^ 2/3[l - (4 H 


^2/3)] +'^32:3, 


(53) 


z' = [aixi + 02X2 + 032^3] 


(54) 


u' = A cos [02;]. 




(55) 



Note that the equations (|48ll53p define three pairs of limit-cycle two-dimensional dynamical systems. The exact 
solutions to these limit-cycle systems are well-known. 
If we define 

Q{t) = aixi{t) + a2X2{t) + asXsit), (56) 

the function Q{t) will be a quasiperiodic time- series. 

Equation (|54[) will provide us with the appropriate nonperiodic truncated exponential behavior. Finally, the solution 
to equation ([55]) will have properties equivalent to these of function ([36]) . The solution to Eq. ([55|) is chaotic in the sense 
discussed in section (|IXp . The maximal Lyapunov exponent of this dynamical system is positive and the solutions 
possess sensitive dependence on initial conditions. 

If we need more variability in the solutions, we can construct a dynamical system such that the right- hand parts 
of the equations depend on function u{t) which is known to be highly nonperiodic. 

These ideas lead to the next autonomous dynamical system: 



x'l = — {xl + 


yl)] -^12/1 + 


eiu, 


(57) 


y'l = 2/i[i - (4 + 


yl)] +UJ1X1, 




(58) 


x'2 = X2[l — {xl + 


■yl)] -^^22/2 + 


£2U, 


(59) 


2/2 = 2/2[l - {xl + 


2/2)] +^2x2, 




(60) 


4 = a;3[i - (4 + 


■yl)] -^^32/3 + 




(61) 


2/3 = 2/3[l - {xj + 


2/3)] +^^32:3, 




(62) 


z' = [aiXi + a2X'2 


i -1- 03X3 -1- a^u 


1 


(63) 


u' — Ax cos[9iz\ -\ 


- A2 cos[92z\. 




(64) 



Now all the components of the solutions to system (|57II64[) are chaotic. 

The solution to Eq. (|64| is chaotic in the sense that the maximal Lyapunov exponent is positive (see section pxp ). 
In the dynamical system (j57ll64|) . the limit cycle subsystems (|57|58p . (|59|60p and (|61|62p are now coupled to the 
chaotic component u{t). 

We now address the function (|37p . We have to construct a dynamical system with solutions that behave like this 
function. 
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Consider the autonomous dynamical system: 





f 2/1)] - 


(65) 


y'l = yi [1 - {xl - 


h t/i)] + WiXi, 


(66) 


x'2 ^X2[l- {xl - 


f2/2)] - W22/2, 


(67) 


Vi = 2/2[l - {x\ - 


1-2/2)] + '^2a;2, 


(68) 


x'z =a;3[l- {x\- 


1-2/3)] -^32/3, 


(69) 


2/3 = 2/3[l - {xl - 


1-2/3)] + '^32:3, 


(70) 


2^4 — 2^4 [1 {x^^ " 


1- 2/4)] - ^42/4, 


(71) 


2/4 = 2/4[l - {xi - 


h 1/4)] + cj4a;4, 


(72) 


- {aiujiyi H 


- 02^^22/2) 2:2, 


(73) 


4 = - {aiuJiyi H 


- 02^22/2) 2:1, 


(74) 


4 = ^ (03^^32/3 H 


- a4tJ42/4) Z4, 


(75) 


4 = - (03^32/3 + aiUJiVi) zz, 


(76) 


u' = Acos[i?iZi + _B4Z4]. 


(77) 



The explanation of this system is similar to that of equations (|48II55I 
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XI. UNIVERSAL FUNCTIONS AND DYNAMICAL SYSTEMS 

Based on the properties of functions pS]). (fTS]). ifTO]). pS]). and ([37|). we propose the following set of equations as a 
universal dynamical system: 



x[ = a;i[l - {xl 4 




(78) 


y[^yi[l-ixl + 


yl)] +cjia;i, 


(79) 


X2 ^X2[l- {xl -+ 


vl)] -<^2V2, 


(80) 


y'2 = Z/2[l - (2:2 + 


yl)] -^1^2X2, 


(81) 


2^3 = 2;3[1 - (0:3 4 


vl)] -^sVs, 


(82) 


y'3 = 2/3[l - (-^3 + 


yl)] + UJ3X3, 


(83) 


X4 = Xi[l - {xl -+ 


-yl)] -t^iVi, 


(84) 


y'i = 2/4[i - (2:4 + 


yl)] +UJ4X4, 


(85) 


z[ = - {aiUJiyi + 


a2^2V2) Z2, 


(86) 


Z2 = - {aiuJiyi + 


a2^2V2) Zl, 


(87) 


4 = - (03^3^3 + 


aiUJiHi) Z4, 


(88) 


4 = - (a3t^3y3 + 


a4a;4y4) ^3, 


(89) 


u = LOinx-y COS 


[-B1Z1 + ^424 + a] . 


(90) 



The solution to Eq. ((90|) is a function with all the properties of Boshernitzan's family of functions (fT9| . Thus this 
family of functions is dense in C[/|. 

The dynamical system (|78ll90p has been constructed using the same technique developed in the papers [1, [l^ about 
universal differential equations. That is, a family of functions known to be dense in C[T\ is utilized as the starting 
point for an inverse-problem procedure that consists in reconstructing differential equations the solutions of which are 
the functions that belong to the universal family of functions. 

This dynamical system can be realized in practice using nonlinear circuits as discussed in Ref . [2lj| . 



XII. CONCLUSIONS 



We have discussed the concept of "Universal functions" and their relevance to the theory of universal differential 
equations. 

We believe that the method of construction of universal differential equations using universal functions is more 
powerful that the method based on splines. 

We have found a connection between the universal families of functions proposed in a very important paper by 
Boshernitzan 13i] and recently obtained exact solutions to chaotic systems. 

We have shown that some functions x{t) that are exact solutions to chaotic systems possess the property that 
yi = x{t) and ?/2 = x{t + r) are statistically independent functions in the sense of M. Kac. 

We have constructed algebraic differential equations that possess solutions with these properties. These equations 
have, in fact, solutions that behave like noise. 

Some known universal equations can only approximate these functions using "solutions" constructed with polyno- 
mial or trigonometric splines. The actual exact solutions to the differential equations are not "noisy". 

We have constructed a system of differential equations, the solutions of which are Boshernitzan's functions. Bosher- 
nitzan's functions are real analytic on M = (—00, 00). One of the families of Boshernitzan's functions are real- analytic 
entire functions on M = (—00, 00). 

We have developed universal- like functions that behave as 5- correlated noise. 

We have constructed physically realizable dynamical systems that possess solutions that are universal-like functions. 

The theory of universal differential equations has been linked from the beginning with applications in analog 
computingd, 0, [nl, [H, [H, [H [3 . 

We believe that the present results can be of interest in the construction of real analog computers because the 
discussed dynamical systems can be realized in practice using nonlinear circuits [21I Isol. Iss!] . 
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They can also find applications in chaos- based secure communications technologies [2l|, [27 
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